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Abstract. In the last months several ballon and 
satellite experiments improved significantly our 
knowledge of cosmic ray (CR) spectra at high energy. 
In particular CREAM allowed to measure B/C, C/O 
and N/O ratios up to 1 TeV/n and PAMELA the p/p 
ratio up to 100 GeV with unprecedented accuracy. 
These measurements offer a valuable probe of CR 
propagation properties. We performed a statistical 
analysis to test the compatibility of these results, 
as well as other most significant experimental data, 
with the predictions of a new numerical CR diffusion 
package (DRAGON). We found that above 1 GeV/n 
all data are consistent with a plain diffusion scenario 
and point to well defined ranges for the normalization 
and energy dependence of the diffusion coefficient. 

Keywords: cosmic rays, propagation model, statis- 
tical analysis. 

I. Introduction 

The problems of origin and propagation of Cosmic 
Rays (CR) in the Galaxy are long standing questions 
which need the combination of several different obser- 
vations in a wide energy range to be answered ffl . 

The most realistic description of CR propagation is 
given by diffusion models. Two main approaches have 
been developed so far: (semi-)analytical diffusion mod- 
els (see e.g. {2| and ref.s therein), which solve the CR 
transport equation by assuming simplified distributions 
for the sources and the interstellar gas, and fully numer- 
ical diffusion models. Well known realizations of those 
two approaches are respectively the two-zone model 
El and the GALPROP 0, @, Q and DRAGON JD 
codes. Generally these models involve a large number of 
parameters which need to be fixed against the observa- 
tions. Their knowledge is crucial not only for CR physics 
but also to be able to constrain/determine the properties 
of dark matter from indirect measurements. However, in 
spite of the strong efforts made on both observational 
and theoretical sides most of those parameters are still 
poorly known. One of the reasons lies in the fact that 
best quality data on CR spectra were available mainly 
at low energy ( < 10 GeV/n). At these energies several 
competing physical processes (e.g. solar modulation, 
convection, re-acceleration) are expected to affect signif- 
icantly the CR spectra to an a priori unknown amount. 



Furthermore, at those energies uncertainties in spallation 
cross section determinations are still sizeable. At higher 
energies, however, only spatial diffusion and spallation 
losses are expected to shape the CR spectra, the latter 
effect becoming less relevant with increasing energy. 
Hence, the study of high energy CR spectra could allow 
to constrain the properties of the diffusion coefficient of 
CR in the Galaxy. 

The scarcity of observational data has precluded this 
possibility for long time. The situation improved recently 
as the CREAM balloon experiment [9| measured the rel- 
ative abundances of elements from boron to oxygen, and 
especially the boron to carbon ratio (B/C), up to energies 
around 1 TcV/n. Furthermore, valuable complementary 
data were recently provided by the PAMELA satellite 
experiment [ 1 1 which measured the antiproton/proton 
(p/p) ratio up to 100 GeV with unprecedented accuracy. 
Also antiprotons are expected to be produced by the 
spallation of primary CRs (mainly protons and Helium 
nuclei). They provide, therefore, an independent test of 
the validity of CR propagation models ifTTl . lfl2l and, 
once these are validated, a valuable probe of dark matter 
models (see e.g. Ifl3l ). 

In this contribution we show that the data recently 
released by the CREAM and the PAMELA experiments 
can fit into a unique plain diffusion (PD) CR propagation 
model. We report the main results of a statistical analysis 
aimed at constraining the normalization and energy 
dependence of the diffusion coefficient of CRs in the 
Galaxy. In order to check the possible dependence on 
low energy effects of these constraints, we also study 
as they change by varying the minimum energy E m i n 
above which data are considered. 

II. Model 

In order to interpret the experimental data we need to 
adopt a theoretical model describing the propagation of 
CR nuclei in the Galaxy between 1 and 10 3 GeV/n. At 
these energies, the propagation of stable CRs is known 
to obey the transport equation |[T4l 

^1 = -V-(D- VN l )+Q t (E k )+ 
at 

- c/3n gas a in (E k )Ni + ^ cf3 'ft" gas 

CTjiNj , (1) 
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where E% = (E — m,A)/A (E is the total energy of a 
nucleus with mass itia — A x m pr ) is the kinetic energy 
per nucleon, constant during propagation as practically 
conserved in fragmentation reactions, (3 is the velocity 
of the nucleus in units of the speed of light c, <jj is 
the total inelastic cross section onto the ISM gas with 
density n gas (r, z) and <Xy is the production cross-section 
of a nuclear species j by the fragmentation of the i-th 
one. We start the spallation routine from A = 64. We 
disregard continuos energy losses, re-acceleration and 
convection, but we check a posteriori the validity of 
this approximation against the experimental data. 

We solve Eq. ([TJ in the stationary limit dNi/dt = 
adopting our numerical code DRAGON |8|. DRAGON 
was validated against well known public codes 
(GALPROP) and against experimental data of sec- 
ondary/primary ratios, as well as 7-ray data. 

We recall below the main assumptions we make. 

A. Spatial diffusion 

We assume cylindrical symmetry and that the reg- 
ular magnetic field is azimuthally oriented (Bo = 
B^lr, z) 4>). Under these conditions CR diffusion out 
of the Galaxy takes place only perpendicularly to Bo- 
Therefore D represents in fact the perpendicular dif- 
fusion coefficient D±. The dependence of D on the 
particle rigidity p is (see e.g. lfT31 ) 

D(p,r,z)=D ex P {|z|/z t } . (2) 

B. CR sources 

For the source term we assume the general form 

Qi(E k ,r,z) = f s (r,z)q , (3) 

imposing Qi(E k ,r Q ,z Q ) = 1. 

We assume that the CR source spatial distribution 
fs(r,z) trace that of Galactic supernova remnants 
(SNRs) as modeled in |[T6l on the basis of pulsar and 
progenitor star surveys ifTTl . 

The injection abundances q Q are tuned so that the 
propagated spectra of primary and secondary species 
match the observed ones. 

For each value of i5 in Eq. (O on is fixed by the 
requirement that at high energy Ek 100 GeV/n, 
at which spallation processes are almost irrelevant, the 
equality a + 6 = 2.7 is satisfied^. 

III. Analysis and results 

Since our main goal is to understand in particular 
the diffusion properties, we want derive constraints on 
Do, 6 and z t in Eq. (0. We consider these observables: 
N/O, C/O, B/C and p/p ratios. They are primary /primary 

'in this regime, the theoretical expectation for the observed flux $ 
on Earth is fa Q{E)/D(E) ~ E- ( . a+s '> 0- 



and secondary/primary ratios. The spectrum of sec- 
ondary/primary ratios allows us to infer information 
directly on S [2 |, but not separately on Do and z t (these 
observables are sensitive to the ratio Do/z t (8)). A way 
to break this degeneracy is to consider unstable to stable 
ratios (e.g. 10 Be/ 9 Be), which are known to probe the 
vertical height of the Galaxy Q. In agreement with lf8~l , 
ifTSl we infer that z t should lie between 3 and 5 kpc. We 
account for a solar modulation potential $ = 550 MeV 
in the "force-free" approximation [19|. 

A. Strategy 

1) B/C ratio: Once the spatial distributions of the 
CR sources and the ISM gas have been chosen, the 
main parameters determining the B/C in a PD model 
are the C/O and N/O injection ratios and the quantities 
5 and D^/zt (which will be always expressed in units 
of 10 28 cm 2 s _1 kpc -1 in this work) in Eq. d2jV 

We fix the source abundances of the oxygen and of 
primaries heavier than oxygen by requiring that they 
match the observed abundances in CRs at E ~ 1 — 
10 GeV/n, while we use primary/primary ratios to fix 
the C/O and N/C0 injection ratios. 

As in [8 1, we accomplish this by sampling, for each 
pair (Do/z t , S), the parameter space (C/O, N/O) and 
computing the \ 2 of our predictions for the C/O and 
N/O modulated ratios against experimental data over the 
energy range of our interest. For the set of parameters 
that minimizes this \ 2 , we then compute the x 2 (which 
we call Xb/c) or " our predictions for the B/C ratio against 
data. By iterating this procedure for several values of the 
pair (Do/z t , 5), we sample the whole parameter space 
of our interest. Minimization of Xb/c l ea ds to m e best 
fit values for (Do/zt,S) and the appropriate confidence 
regions. 

2) Antiprotons: The construction of a statistically 
meaningful variable for the p/p ratio is rather simpler 
than for the B/C. Indeed, if we neglect the systematic 
uncertainties associated to the production and interaction 
cross sections, the propagation of secondary antiprotons 
depends essentially on Dq/zu S and the source abun- 
dance ratio He/p. This last unknown can be easily fixed 
by looking at the measured spectrum of He at Earth, 
which is relatively well known. Therefore, we construct 
a Xp/p by comparing our predicted p/p spectrum for dif- 
ferent values of (Do/z t , 8) to experimental observations. 

3) Joint comparison: Since we have two different 
data-sets for the same physical framework, it is impor- 
tant to find links for joint optmization of the combined 
analysis. The two data-sets being uncorrelated, it is 
possible to define a joint x 2 knowing the single results 
of the previous analysis 

X 2 = \ (x 2 b/c + X 2 p/ P ) (4) 
and minimize it with respect to (Do/zt, 5). 

2 Note that N = 14 N + 15 N is a combination of primary and 
secondary nuclides. 
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The complementary nature of the two data sets will be 
demonstrated in section IIII-CI in which the overlapping 
of the separated analysis is evident. 

B. Experimental Data 

So far the best B/C measurements above 1 GeV/n 
have been provided by the HEAO-3 |2__ and CRN ED 
experiments in the range 1 < Ek < 30 GcV/n and 
70 GeV/n < E k < 1.1 TeV/n. Recently, the CREAM 
[9 1 experiment has released data l22l improving signif- 
icantly the available statistics at high energy. C/O and 
N/O data are taken from the same experiments as well. 

For antiprotons we use experimental data released by 
BESS for the periods 1995-97 |2__ and 1998 El in the 
energy interval 1-4 GeV, and by CAPRICE (1998) |2__ 
in the range 3-49 GeV. Recently also the PAMELA 
experiment has released p/p data in the energy range 
1 — 100 GeV ITOl . We include them in our analysis. 

Solar modulation has been demonstrated to have im- 
portant effects in the determination of the p/p ratio at 
low energy. We account for modulation in the force 
free approximation. For each data set we use the solar 
potential of the year when data were taken. 

C. Results 

In Fig. _] we show our results for the separated and 
joint analysis of the B/C and p/p data. The left-hand 
plots represent the 1, 2 and 3a contour levels of the 
Xb/c ™ ^ e (Do/ z t, S) space, the central plots show the 
same contour levels for Xp/ P > while the right-hand side 
plot is the result of the joint analysis described in section 
ITILA3? . From the panels in the first column the impact 
of CREAM results on our knowledge of the propagation 
parameters is evident: CREAM high energy data favour 
a smaller value of S with respect to previous ones. 
Also PAMELA p/p data are sensitive to our propagation 
parameters, in contrast with pre-PAMELA p/p data (see 
the upper row panels in Fig. [TJ. Finally, it is clear from 
the lower row panels as for E m i n > 5 GeV there is an 
almost complete concordance between the CL regions 
constrained by high energy B/C and p/p data. Indeed, 
even the la regions do have a significant overlapping 
and the minimum joint \ 2 is 1.2 for Bmin = 6 GeV/n. 
The joint analysis indicates a best-fit value of 5 < 0.5, 
possibly favouring a Kraichnan power spectrum for the 
turbulent galactic magnetic field. While the best fit for 
S seems not to be strongly dependent on E m i n , we 
find that the best value of D$/zt tends to be larger 
when a larger .Emin is considered. This could indicate, 
in agreement with naive expectations, that the scale of 
vertical diffusion z t is smaller at higher energies. 

The concordance between nuclear and antiproton data 
is also evident from Fig. [2] where the predictions of the 
combined analysis best-fit model are compared with B/C 
and p/p experimental data. 

Finally, in Tab. _] we recap the findings of our min- 
imization strategy. Best fit values for (Do/z t ,5) are 



obtained considering first the two data sets separetely 
and then jointly. 

Our analysis shows clearly that the two data sets are 
statistically compatible, within our model. 

IV. Conclusions 

We performed a statistical analysis to constrain the 
CR propagation parameters in a plain diffusion sce- 
nario numerically implemented in DRAGON. Taking 
advantage of the new CREAM and PAMELA high- 
energy data we performed a combined analysis of B/C 
and p/p data in several energy ranges. This approach 
allowed us to test the (weak) dependence of our results 
on poorly know low-energy physics. We showed that 
above few GeV/n the whole data sets considered here 
are consistently reproduced by a PD model. Our findings 
favour a Kraichnan type (S ~ 0.5) dependence of the 
diffusion coefficient on energy. 
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TABLE I 

Results of the statistical analysis described in Sec.HTTI 
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Fig. 2. Comparison of experimental (top of the atmosphere) data for B/C (left panel) and p/p ratios with our joint analysis best-fit model 
(second row in Tab. U. Dotted lines refer to LIS ratios, accounting for solar modulation with potential <j> = 550 MV. 



